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ABSTRACT 

A formula is derived for the mutual inductance of two helical wires wound 
bifilarly upon the same cylindrical form. This aims at a precision sufficient for 
use in absolute electrical measurements, so that terms are retained which, com- 
pared with the principal term, are relatively small quantities of the order of the 
square of the ratio of pitch to cylindrical radius. Nominally the two helical 
windings consist of wires of equal diameters, equally spaced, with equal cylindrical 
radii, but allowance is made for the case where these three equalities are not 
quite attained. The axial lengths are assumed to be equal, so that, since their 
pitch must be the same, the total number of turns is the same. The latter is 
considered to be integral, which permits of the two return currents being carried 
each by a straight wire parallel to a generator of the cylindrical form and close 
to it. The effect of these lead wires is included in the formula. The azimuth 
angle between these generators is assumed to lie between — w and 7r, so that the 
point of current entrance for the primary is on a turn adjacent to that of the 
secondary. Allowance is made for a general type of current distribution. 
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I. INTRODUCTION 

Mutual inductors are now being constructed as single-layer, bifilar 
windings on a circular, cylindrical form. The dimensions may be 
measured as precisely as those of a self inductor. However, a for- 
mula for computing the mutual inductance between bifilar windings 
is not available. The object of this paper is to supply that need. 
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II. MUTUAL INDUCTANCE, M, OF THE WIRES AS 
AN INTEGRAL OVER THEIR SECTIONS 

The helical wires, w x and w 2} have the pitch 2irp. The two full cir- 
cles of figure 1 with diameters d x and d 2 are their generating circles, 
for these circles sweep through the volume of the wires when the 
plane of the paper is given an angular velocity, ddjdt, about the x axis, 
together with a translational velocity, pdOjdt, in the ^-direction. The 
same definition applies to a helical wire with any generating section, 
such as an ellipse or rectangle. The modifications necessary to extend 
the formulas to wires with general types of cross section will be indi- 








x=x°-2irp 




X«TTp 



M. 



x«zflp 



Figure 1. — Reference circles which generate the two helical wires. 

cated later. The center of the first circle is at z— 0, y=a u 2=0, that 
of the second at x=z% 2 , y—a 2) 2=0. In the actual construction, the 
attempt will be made to have the two wires equally spaced (x% 2 =Trp) 
and to give them the same mean cylindrical radii (a l =a 2 ). As neither 
of these conditions may be attainable with sufficient precision, we con- 
sider the more general x% 2 , and assume a x —a 2 small compared to the 
pitch. The two adjacent circular sections bounded by full circles in 
figure 1 may be designated as the "reference sections" of the wires, to 
distinguish them from the initial or end sections, where the current 
enters or leaves the wires. However, the first of these sections, Si, 
may be considered the initial section, 8 %v of wire w u but the second, 
S 2 , is in general not the initial section, S iv of wire w%. The initial 
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section, S i2 , is the circular section in which w 2 is cut by a plane which 
passes through the x axis and makes a given angle, a, with the plane 
XOY, where — t ^ a < it. 

If x, y, z are rectangular coordinates and x, r, the cylindrical coor- 
dinates with the axis, as shown in figure 1, then the azimuth of the 
initial section, S iv of w x is 0=0^=0, and that of the end section, 
S ev is 6=d ei =2irN, where N is a (large) given positive integer, the 
total number of complete turns. The initial and end sections thus 
being in the same plane, the lead wires for the return current may be 
made to consist mainly of a straight wire parallel to a generator of 
the cylindrical form and close to it. Similarly, for wire w 2 the initial 
azimuth is 0=0^=^ and the end azimuth is 6 e2 =a-\-2irN; so a simi- 
lar lead wire in the plane d—a. may be used. These are shown, to- 
gether with the central helical filaments of W\ and w 2 , in figure 2, 
which is drawn for the case where the wires are equally spaced (x% 2 =Trp), 
and where the azimuthal angle is a= —w (as explained below). 

To designate any particular helix, h x , of the helical wire, w u it suf- 
fices to assign the two coordinates, x\ and r u of the point P? in the 
reference plane of figure 1, through which this helix passes. Then, if 
%u Vu z i ar e the rectangular coordinates of any point P x on this helix, 
the equations of the helix, h Xy are 

Xi=x° 1 +pd 1 ] 

y x =r x cos 6A, where 0^d l ^2wN. (1) 

Z\=Ti sin Q x \ 

It is understood that x\ and r x are constant (as well as p), and that 
as the line parameter, $ u ranges over the prescribed interval, the 
point P x moves along the helix. The pair of constants (x° ly r x ) may 
have any values subject to the limitation 



af+Cn 



-KfJ 



The particular pair, x?=0, r x =a x , corresponds to the central helix of 
wire w x . Its equations are 



x x =pd x 
y x =a x cos 0i 
z x =a x sin 2 



, where 0^0^ 2^. 



Similarly, any helix, h 2y of wire w 2 , may be specified by assigning the 
coordinates of x? 2 and r 2 of the point P% in the reference plane of figure 1, 
through which this helix passes or would pass if continued outside 
the range which limits the line parameter, 2 . The equations of h 2 are 



X 2 =xl+pe 2 =x i2 +p{0 2 —a) 
y 2 =r 2 cos 2 
z 2 —r 2 sin 2 



, where <d 2 <a+2irN. (2) 



This helix passes through the point in the initial section, S iv whose 
cylindrical coordinates are (x i2 ,r 2 , a), where x^x^—pa. The pair 
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of constants (xl, r 2 ) which characterize h 2 are subject to the inequality 



(x° 2 -x° 2 y+(r 2 -a 2 y 
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The particular pair, xl=xl 2 and r 2 =a 2 , correspond to the central helix 
whose equations are 



x 2 =x o 2 +p6 2 
y 2 —a 2 cos 2 
z 2 =a 2 sin 6 2 



, where a^d 2 ^<a+TrN. 
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In the special case where the helices are equally spaced, (x$ 2 — Tp), 
the central helix passes through the point in its initial section, S i2y 
whose cylindrical coordinates are (p {-*-{- a), a 2 , a) which, for the case 
a= — t, becomes (0, a 2 , —w). 

This is the point 2 ' in figure 2. The two wires then have the same 
mean end planes. In the general case, if — 7r<a<0, the point O2 
would be moved toward 2 along this helix. If 0<«<V, it would fall 
between 2 and 0". 

Every helix in either wire has the same axial length, l—2irpN. 
Every helix, hi, in the wire w x has the same initial azimuth, 0^=0, 
and end azimuth, 6 ei =2irN. They differ in their cylindrical radii, r u 
and in the initial and end planes. Thus the helix hi lies between the 
initial plane, x=x h =Xi, and the end plane, x=x ei =Xi+l, these limit- 
ing planes being and I for the central helix. 

Every helix, h 2 , of wire w 2 has the initial azimuth, 6 i2 =a, and the 
end azimuth, 6 e2 =a-\-2wN. The helix h 2 lies between the initial 
plane, x=xi 2 =x 2 ) -\-pa, and the end plane, x=x e2 =Xi 2 +l, these limit- 
ing planes for the central filament being x° Q2 and x% 2 -\- 1. 

It is assumed that I and a x (or a 2 ) are finite quantities of the same 
order of magnitude, and that the pitch, 2 7^, is a small quantity of the 
first order compared with a x or a 2 . A formula for the mutual induc- 
tance of the two wires will be obtained which neglects terms that are 
smaller, relative to the principal part, than the square of this infini- 
tesimal compared with one. It is evident that the diameters of the 
wires must be small quantities of the same order as the pitch. 

In the formula to be obtained, the pri ncip al term is proportional to 
di^lv 2 } s0 t na t a term proportional to ^a x a 2 fp would be (relatively) a 
first-order small quantity, the proportionality factors being "finite" 
functions. Hence a finite function would be relatively of second order. 
Such terms are not quite negligible, but they need not be evaluated 
with great precision. In fact the principal term must be of the order 
of magnitude of 10 mh (10 7 cm) or larger, if the coil is to be of value in 
absolute electrical measurements, so that in the "finite", or smallest 
terms, retained in the formula, an error of a/2 is tolerable. The 
mutual inductance of a bifilar coil is practically equal to the self 
inductance of one of the windings. In the case of a coil of 500 turns, 
100 cm long with radius 14 cm, the total self inductance is 17 (10) 6 cm, 
so that a/2 is 7 cm, or about 4 parts in 10 million. For a shorter coil, 
27 cm long, 135 turns with a larger radius, 22 cm, the self inductance 
is 7 (10) 6 cm, so that a/2 = ll cm is about 4 parts in 3 million. This is 
the basis for replacing certain expressions defined as infinite series or 
as definite integrals by simpler expressions in closed form, which are 
sufficiently close approximations. 

Concerning the distribution of the steady current in the wires, we 
make the assumption that the vector-current density at any point, 
Pi, of wire w% has the direction of helix h x passing through that point, 
and its magnitude remains constant as the point moves along this 
helix, so that it is some function, Ui(tf, ri)=u(xi—pd u ri). The 
current density will then be a solenoidal vector. Similar assumptions 
as to the current density in w 2 make its magnitude some function, 
u 2 (x% r 2 ). 
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Any axial plane (0= constant) cuts the wire w t in a circular section 
of diameter d x . If dS x is an element of area of this section, a helical 
tube of current, whose (oblique) section is dS x , carries the current 

since [l-\-p 2 /ri]~* is the cosine of the angle between the normal to 
dSi and the direction of flow. The total current carried by the wire 
is unity, so that 

* u x (x° u r x )dS x _ 



ip- 



where the integration may be taken over any section of the wire made 
by the plane 6= constant. 

A precisely similar equation applies to the unit current in wire w 2 . 
It will be convenient in these integrals, and in those over the same 
sections which are to be encountered presently, to consider them as 
taken over the circular reference sections in the plane of figure 1, 
even though in the case of the section S 2} this may not be an actual 
section of the wire w 2 . 

If the wires were of rectangular section, it would be appropriate 
to take dS x =dx° x dr x and dS 2 =dx 2 dr, but since they are circular, it is 
more natural to use plane polar coordinates (p u <t> x ) for points in 8\ 
and (p 2 , fa) for points in £2, where 



while 
so that 
where 



x\=pi cos 0i and r x — a x =p x sin <j> x 
A~ x o 2 — P2 cos <f> 2 and r 2 —a 2 =p 2 sin <j) 2 , 
dSi = p\dp\d<t>i and dS 2 = p 2 dp 2 d<j) 2 j 



If rn a (x u r x ; x%, r 2 ) denotes the mutual inductance (Neumann's line 
integral) of the two helices, h x and h 2 , then the mutual inductance, 
M a , of the two helical wires, w x and w 2) is 

,, f f tti(s?, r x )dS x f f majx !, n; x° 2 r 2 )u 2 {x° 2y r 2 ) ^ f . 

Ma= i J ir+mi J vt+f^ **■ (4) 

We shall obtain an expression eq 52 for m«, whose principal part is 
constant; so this term contributes its own value to M a when used in 
eq 4, as shown by eq 3 for any shape of wire section. The smallest 
terms of m (second order) may then be evaluated by taking 

Ux _l^ , _ u 2 1 
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where these terms are concerned in the integration of eq 4. This is 
also valid for any shape of the sections. For the small terms of the 
first order we use (when the sections are circular) the Taylor expansions 

75 1/ o =-q h+fr|Pi cos 4>i + c n Pi sin 0! 

u (x° r ) 1 f I 

/! 9/ , = gr U+^ 2 P2 cos <fe+c r2 p 2 sin & > 

where the constants c H and c ri are the values, at the central filament of 

c *i=^D x oUi and c Tl =-D r Ui (5b) 

with similar expressions for c X2 and cv 2 . In the case of most practical 
importance, the two wires will be of the same diameter and will have 
received the same treatment, so that c xl =c X2 =-c x and c n = c T2 =c T . 
In this case it will be found that the variation with x° is unimportant, 
as the effect in one wire cancels the effect in the other, so that c x need 
not be known. The radial variation, however, has some terms which 
cancel and some which add, giving a term of second order involving c r . 
It may be noted that c r is zero if the current is uniformly distributed 
over the cross section of the wire. In the "natural" distribution 
where the current density varies inversely as the cylindrical radius 
(that is inversely as the length of helical filament), the value of c T is 
— 1/a, where a is the radius of the central filament. If the con- 
ductivity of the wire were uniform, it can be shown that the distribu- 
tion would be "natural." 

It may be noted that the integrations in eq 4 are taken over the 
plane reference sections of figure 1. The equation is valid for wires 
with the most general type of cross section. If the sections are not 
circular, the second members of the two eq 5a must be replaced by 
other first-order approximations with other variables of integration 
which are more suitable for the particular type of sections considered. 
These equations are only necessary for use with the small correction 
terms which are of first order. 

We leave eq 4 at this point and, in section III, digress with tedious 
detail in pursuit of an approximation for m which may be used in 
the eq 4. The results are collected in eq 53. After that, in section 
IV, it is relatively easy (with circular sections) to evaluate the funda- 
mental integral given in eq 4. 

III. MUTUAL INDUCTANCE, m, OF TWO ELEMENTARY 
HELICAL FILAMENTS 

1. FORMULA FOR zn IN TERMS OF THE co-FUNCTIONS 

The general formula for m has been developed in a recent paper, 1 
in which the small correction terms depending upon the azimuthal 
angle were reduced to a practicable form by limiting the application 
to cases where the two helices were nowhere close to each other. 

* C. Snow, Mutual inductance and force between coaxial helical wires. J. Kesearch NBS 22, 239 (1930) 
RP1178. 
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These terms require a different treatment in the present problem of 
the bifilar inductor, but for all other terms we merely quote the 
results of that paper. For brevity, this will be referred to as paper 1. 
It is first necessary to place here for reference the function of the 
modulus k denoted by <t> n (k), which is in fact equal to 7r -1 Qn-\ with 
argument (2— & 2 )/& 2 , where Q is Legendre's function of the second 
kind. Since 0<fc<l, it is defined for n=0, 1, 2, 3 * * * by 
the hyp ergeome trie function 

^ k) = k Zr(2n+l) F{n+h n+h to+U.P) (6) 

or in terms of the complementary modulus £'=(1— -A: 2 )*. 

Wfl+lf 2 

*.(*)=- V"P (7i+ *' n+h 1; k) log v 

+ f ^p^ T -^^[Hs+n + i)-Hs+im t (7) 

where \f/(z) is T'(z)/T(z). The first two functions are given in terms 
of the complete elliptic integrals 

7rcl> (k)=kK(k) and T4> l {k)^[K{k)-E{k)]-kK{k). (8) 

When k x is a function of x, r u and r 2 given by 



then 



_4rir L _ J a 2 +(r 2 -r 1 ) 2 

* x 2 +(r 1 +r 2 ) 2 ° r ^ 4nr 2 +x 2 + (r 2 -n) 2 ' w 

00 

<t>n(k x ) = ^Wjf-W'Mr^Jn^ds, (10a) 

where J n is BessePs function. 

Also when k is given by eq 9, <f> n {k) is a solution of the partial 
differential equation 

D\+D\ + 1 ^D ri +^^^H n (k)^0 ) (10b) 

where a is any constant. The same equation holds when r x is replaced 
byr 2 . 

The function <t> n (k) is the coefficient in the following Fourier's series: 

00 

^l-t>eos>e^ Mk)+ 2j» ik) C0S 2n6 - 



n=l 



In the following, it will be assumed that k x and k' x are the functions 
of x, ri and r 2 defined in eq 9. When x is merely replaced by x lf this 
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is indicated by k\ ; when x is Z, by k t ; and when x is zero, by k Q ; similarly, 
fexj when x=X. The fundamental equation is eq 13 of paper 1. 

m=o)(x e2 —x il )—o)(x e2 —x ei ) + u(x i2 —x ei )— o>{x i2 —Xi^ (11) 

where the subscripts i and e indicate initial and end planes, and co(x) is 
written in place of a)(x, r x , r 2 ; a). This function was there resolved 
into three constituents, as follows: 

co(x)=W6(x) + a) x (x) + QOa(x,<x). (12) 

The term co a was called the azimuthal co-function, as it is the only one 
depending upon a. All three functions depend upon r x and r 2 . 

The associated current sheet of the helix, h u is defined as the circular 
cylindrical surface coaxial with it, having the same radius, r i} and end 
planes, x t , x e , on which surface the linear density of current has the 
angular component, jd = l/2Trp=n=th.e number of turns per centi- 
meter. The axial component of current isj x =l/2Tvr u so that the unit 
current in the x direction is uniformly distributed over the sheet. 
The two terms, o)e(x)+o) x (x)j when used in eq 11 give the mutual 
inductance of the two current sheets which are associated with the 
helical filaments, hi and h 2 . The first, we(x) } is that due to their 
angular or 6 components of current, and this contributes the principal 
part to m. The second, u x (x), is relatively a second-order inflnitestimal, 
representing the effect of the axial components of current. Each of 
these functions is an even function of its z-argument, vanishing with 
it. This property used in eq 11 gives the following interpretation 
of these functions. 

Let the helices have the same end planes, x { =x t and x e = Xe . 
Then, if x=x t — x e , the function 2ooe(x)-\-2a) x (x), or explicitly 
2ooe(x, n, r 2 )+2co z (x, r ly r 2 ), represents the mutual inductance of two 
coaxial circular cylindrical current sheets of axial length x and radii n 
and r 2 , when they have the same end planes. 

The integral definitions of these functions are 



dx 2 faiki 
Jo Jo 



I f*x rxt 

«*(^)="7= = d** I #o(&i)(fci ( 13d ) 

V^2jo Jo 

The azimuthal function, co a (z, a), is written 

w a fe a)= S «%fo a )> (13c) 

71 = 1 



where 




w n (x, a) = 


_<y/r 1 r 2 

v 2 



rx rx 2 

dx 2 [<t> n - 

J0 JO 



i(*i) + *«+i(*i)] coz-fa—x+pofidxi. (13d) 

p 
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It will be found that a) a is also infinitesimal compared with the 
principal term co§ but, unlike co^, it contains first-order terms. 
If, for brevity, we define the small quantity X by 

\=x i2 — x tl =x%— xl+pa, (14) 

the fundamental eq 1 1 becomes 

m=«*(Z+X) + tt*(Z— \)_ 2a (\)+2a> x (l) 

+ co a (7+A, a) + co a (— Z+X, a)— 2co a (X, a). 

Since X is small, this becomes 

i»=2« tf © — 2«# (X) +2« X (Z) +4Vv 2 (2~) 2 7r 2 1 (*,) 

+ co a (Z+X, a) + co a (-Z+X|a)-2(o a (X, a). (15) 

In the evaluation of these co-functions which follows, frequent use 
will be made of the results and methods of an earlier paper, 2 which 
will be referred to as paper 2. Reference will also be made to paper 
3. 3 

2. PRINCIPAL TERM, «,(J) 

The principal term, 2c*o(l), in eq 15 represents the mutual inductance 
of two coaxial current sheets of radii r x and r 2 , which have common 
end planes and axial length, I. It is due to their angular component 
of current only. When r 2 — >ri, this becomes the self inductance of the 
circular cylindrical sheet. 

Since r 2 — /^ is small, it will be appropriate to express this term 
as a self inductance corresponding to the cylindrical radius a, which 
is the arithmetic mean of the radii a x and a 2 of the two central fila- 
ments. 

The self inductance of the current sheet of length, I, and this 
radius a is 

£.=^VF+^[#-#fT? (#-*)]> (16) 

where the modulus k is now given by 

k 2 =j2^2 and D=2a=a x +a 2 . (17) 

In eq 19 of appendix 3 of paper 2 is given an expansion which, to 
the second order inclusive, is, by a slight change in notation, equiva- 
lent to the following equation: 



1 C. Snow, Formula for the self inductance of a helix made wilh wire of anv section, BS Sci. Pap. 21, 431 
(1926-27) S537. 

3 C. Snow, A simplified precision formula for the inductance of a helix with corrections for the lead-in wires, 
BS J. Research 9, 419 (1932) RP479. 
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where the constants A%, A Z} and A4 are functions of thejmodulus, k, 
which was defined in eq 17, as follows: 



4.-CH 

. _ir ±E-(?-k 2 )K 5~1 
^ _ 8L~ ~T~ 2 J 



(19) 



It is evident that in the fourth term of eq 15 we may place V r i r ?— a 
and (t>i(ki) = <t>i(k) 7 since this term is of the smallest order retained 
(second order). 

3. SMALL TERM, w#(X) 

To simplify the term — 2co^(X) (on account of the fact that X is 
small), we write, by eq 15 and 13a, 



2V0 



> 2 Jo 



-2«i(X) = ^-/- 2 (\-x l )4>i{h)dx l . 

V Jo 

But since X is small, Xi is small throughout the range of this integration. 
Also r 2 — r x is small, so that by eq 9 k[ is small. Therefore, in this 
integral we may use eq 7, which gives 

— 20 1 (fc 1 )=~ log k[+\p(^j—\P(l) + infinitesimals 
= ^[log£;+2-2log2] 

=i{4-2 log 8V^2+ log [xt+ (^-n) 2 ])- 

Hence co*(X) is an even function of X, which differs infinitesimally 
from the following integral: 



-2coa(X)=^ f X (X-*i)[4-2 log Sjw + log K+fo-r,) 2 ! 



dxi. 



The result of this integration and of others which are to follow may be 
abbreviated by introducing the complex variables 

^HX^i^-^^-^-^+iixl-xl+pa). (20) 

and 

2 > ^\r 2 - ri \+i\\\. (200 

Then ^(X) is a relatively small term of second order, given by the 
real part of 
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4. AXIAL-CURRENT TERM, co x (l) 

This term of second order may be treated as a constant. It may 
be computed by 

2^(0=2^(0-^(1-1), (22) 

where, by eq 43 of paper 1, 



■®=** l± ^+w& 1 °**- 



(23) 



This approximation was derived on the assumption (here adopted) 
that terms in the co-functions of order of magnitude .12? are negligible 
(3 parts in 10 7 ). 

In paper 3 (eq 8) an approximation was used which is equivalent to 

(k'\ h 2 dk f 

^-*-logf ,if*'<*. (23') 

Either eq 23 or 23' is adequate for computing this relatively small 
term. The former has the advantage of a single form, and is one 
which blends naturally with the corresponding terms representing the 
effect of the lead wires (eq 62, below). 

5. AZIMUTH TERM, w a 

In paper 1, pages 255-258, it is shown that instead of the integral 
13d for u n we may take 



[<01 

Tnx ( t\ ', /j ; n\r a -n\ \ . / x\l 

: — cos ni a )+( — — '— a ) sin nl a ) \e p 

\ C a se-\ x \ s J{n, 
"Jo [W^ 



«N- ri| 



™ Jo r i+ (|A 2 T*L (24) 

where for brevity 

J(n, s) m J n _ x (r x s) J n _! (r 2 s) +J n+1 (ns) J w+1 (r 2 s) . (25) 

The upper or lower of the ambiguous signs is to be taken according as 
x is positive or negative, the quantity affected with the ± sign being 
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zero when x is zero. This expression 24 is not equal to the integral 
definition of o) n given in eq 13d. It is, however, effectively the same 
since terms have been discarded which cancel out of eq 11. Such a 
procedure is allowable if we make eq 24 the basis for simplifying the 
function co n (k) both in the case where x is small and when it is finite, 
so that the terms in question actually are discarded from each of the 
four co-functions occurring in eq 11. 

The simplification of eq 24 depends upon the fact that we may 
neglect all terms in this equation which are infinitesimals of higher 
order than the second in comparison with coa, which contains the factor 

When x is finite, the second integral in eq 24 is negligible, and the 
first may be replaced by 

m 
JO V r l r 2 

as shown by eq 10. This is an even function of x. 
Hence 

u a (l+\ a)-\-a) a (l—\ a) = 

-n\rt-ri\ 

cos na 



, — J2ZVM -% — 1 n{x° 2 -x< {) 2VV ..... n w 

=w^^2j» cos v ^2j [0n " i( ^ +0n+i( g)] 

n=l n=l 

= - Wnr 2 |— log [1 - e { p J] 

+ | ^-i(*i) + ^(t3] cos , a j, (26) 

where the real part is to be taken. 

The logarithmic term containing the factor l/p is the most important 
correction term. The series in cos na is the smallest order of magnitude 
retained (second order relative to &$). It is shown in eq 57 of paper 1 
that an approximation in finite form sufficient in all cases is given by 






-(|) 2 (co S f)sin-(fcco S f)-^ 2 sinflog ^inf+^/l-Fcos^ l (2?) 

(By a misprint the logarithm had the factor \\k instead of 1/k 2 in 
eq 57, but this was used correctly in the equations depending upon it.) 
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The derivation of eq 27 depends upon the fact that the continuous, 
even function of 0, f(d), which is never negative and is a periodic 
function of 6 with period 27r, defined for all real values of by the series 



- /zA 2 4 VV cos nd 
is represented approximately, for all values of 0, by 



M- 



sm 2 



The exact expression for f(d) is 

/W =4 (^:Y 1 -^:) for the ran S e O^^* 



or 



™={\£\-(lJl ioTthe 



range — 2ir<0<2ir. 



(28) 



(280 



(28 ,/ ) 



(28'") 



The graph of /(0) is a continuous curve represented by a succession of 
parabolic arcs. Its slope is discontinuous at 6=0 or 2irn. 
Using eq 27, we find that expression 26 becomes 

u a (l-{-\ OL)-\-a> a (l—\, a) 

t/ lrt-ril-K(x;-i?) \n 
1-6 J 



1 



+T* l( ^ )+ itW 1 "^ cos2 f 

^(IJ'G 08 !) 8111 " 1 ^^ ^) 

a &,sin|+-Wl — /^cos 2 | 
f sin j. log[ \, 



(29) 



We next consider the last term in eq 15, which is — 2w a (X, a). Since X 
is small, we must now consider all the terms in eq 24 in order to obtain 
valid expressions for co n (X, a). Considering X as an infinitesimal of 
the same order as the pitch or smaller, it is evident that the ratio, 
X/#, may be infinitesimal, but in general may have a finite value 
positive or negative. 
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Reference to eq 10a shows that we may write 
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- lb 



-Cs)' 



+(? 



+ 



Vnra 



»)]' 



/*oo 



«-| x I* J(n, s)cfe= 



1 



[«-i(ix)+*H.i(Jfcx)] 



['+(? )"J 1+ (f )■ 



6- ,xl *J(n, *)<&, (30) 



where the last integral is no longer infinitesimal with p, because X is also 
infinitesimal. In fact, the substitution s=nt/p gives 



7=7^2 



e-^'Jn, s)ds 



which shows the integral is a function of the finite ratio, |X/p|. 
Referring to eq 25, we find 

<«< 1)= J --(?<M?0 W?<M?<) 

in which the argument, ntr/p, of the BesseFs function is large (even 
in the case n—l), except when t is small, but in that case the factor t 2 
in the integrand of eq 31 renders its values of little or no importance 
in the initial range of the variable t. (Also all the BessePs functions 
involved vanish with t except J ). We make an infinitesimal error in 

I by introducing in place of jl n, —) the principal term of its asymp- 
totic expansion, which is 

l^j(n, |)*|[cos|(>y-f a )t-(-l)" sin ^M-r,)*]- (32) 

Reference to tabular values of BessePs functions shows that this is 
exceedingly accurate for arguments greater than 10. In this case, 
if the pitch, 2rp, is 0.1 cm, and r 1 =\b cm, the argument nr\t\p is 
3007r?2i, or £(10) 3 , even in the case n=l. Hence the asymptotic 
expansion is good if t is as small as .01, in which case the factor t 2 of 
the integrand in eq 31 is not greater than (10)~ 4 . 



229032—40- 
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Hence, if X, r% — r% s and p are all infinitesimals, the integral / of 
eq 31 may be shown by use of eq 32 to differ by an infinitesimal 
amount from the following: 



^MaW+FR 2 }"^" cos I ^- r ^ tdt 



- -*w, 



This is equivalent to 

/" . inz' 
-[&"■*. < 3 » 

where the real part is to be taken and z' is defined by eq 20'. 
The detailed justification of this process of deriving eq 33 is given in 
appendix 5 of paper 2. In appendix 8 of the same paper, it is also 
shown by deforming the path of the integral that 






Hence eq 30 becomes 



^r x r 2 I 

* Jo fl- 



■(f)" 



XU 2 e~ W ' J(n, «)&=~[^- 1 (Ax) + ^fi(tx)] 

n ./ /\ ««i o ■ tcosnt—n( — ) sin nt 
+ | + V 1+ «£V-^ + 2 I W * (35) 

By the same process it is found that 

se-M'Jjn, s) , 2^ | " p 7, 

fi+^YT ™ 1+ * 2 

V* i nz'\ — Ziz P co smnt— nt cos nt ,. /rtn . 

= ( H )^ » 2~ I 7-7Y2 «*. (36) 
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Now, by eq 13c 
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7=- — 7= / Wn{\ OL), 
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where co» is given by eq 24 and where the two integrals of eq 24 are 
given by eq 35 and 36. Hence 



w fl (X, a) ^cos na 

7r V r l r 2 n = 1 ^ 



~[<t>n-l(h)+ct>n + l(h)]+^ + 



+^+f) e --?+IX' 



/ cos nt- 



-JiJ si 



sin nt 



dt 



'+(1) 

l/\ . W2'\ ~F~ • 2is sin na r°°smnt—nt co 
— -( l-\ )e v smna , I ^-7^ 



njra-nl 



co ^4 



where the real part is to be taken, z' being given by eq 20'. The upper 
or lower of the ambiguous signs belongs with X positive or negative, 
respectively. 

Making some cancelations, this gives 

co a (X, a) __ 2^ COS na 1" (±\i A n. m cos na 

ir-y/rfa * n=i n Tr n =i n 



+ 



^%!'S 



2 cos n< cos na 






2 sin nt cos na 

71 



i^T " 2 sin ^ s i n nQ , • 2 cos r?i sin na"]} , nn , 
± 7LS «? '§ n 1 (37) 

It is now easy to see, by taking the real part of this equation, that 
the ambiguous sign ± may be avoided, for we obtain the same real 
part if we replace z' by z in eq 37 and take the upper of the double 
sign T . We thus return to the complex variable z of eq 20 and have 
no more need for z' . 

It is necessary to evaluate eq 37 for the case where — 27r<a<^27r. 

Writing 2 cos nt cos na=cos n(t+a)-\-cos n{t— a) y and then making 
the substitution t'=t+a and t f =t—a, respectively, we obtain integrals 
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over the infinite range, in which the integrand contains either the 
factor 

CD 

7r~ 2 2^~ 2 c os nt 
i 

or 

CO 

w^y^.n' 1 sin nt } 
1 

both of which are periodic functions of t with period 2ir, the first being 
an even function of t, which is equal to 



\2ttJ 2tt^6 



for the positive range, 0<£<2tt; the second an odd function, which 
is equal to 



-hliih] 



for this range. 
Hence, writing 



J 1 » co /»2x(n+l) 

(...)*=2 (•• 

«=0 J 2m 



)dt, 



and taking account of the periodicity of the series, we find (for the 
case 0^|a|<27r). 



7 * f -_*_LV: 



2 co s nt cos na z 2 \\ 2 sin nt cos na y 

n 2 v 2 / i n 

i J 



, f" Q-IX"-'+l>fe->-» 



V 

Vf |^<)('-4feK " 



n=0 



Snow] 



Bijilar Mutual Inductance 

('+"-&X''-'+l) + fe)' (2, - i) 
(<+*-£M£)' 
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(38) 



and 



t =:=*? C dt IV : 

J. PI l 



n* 



'2- 



2 sin nt sin ?&« /ST^ cos ^ s i n wa 






0-S*fe) 

'f-t+^-((+»+£)(2l-D 

(<+»+£)'+fe)' 



P-| + |-( (+ »-^)(2l-l) 



("HSW^ I 



(39) 



where the upper sign of the T is to be taken when a is positive and 
the lower when a is negative. 
Hence 



li+Ji=-2 



I s 



Mfl+H2JS> Ti > 



*- 



27T^ 



eft 



+ 



Sf 



i* 



* , - < +H^) (? '- l) 



t+n- 



2rp 



fa-\-iz 



(40) 



t+n- 



2irp 
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or 



/i+/a=-2 



n=0 



Jr*a/2x 



'^i+i-^ 



h§+(i)' + ^ 



\2wpJ 



pa+iz 
2irp 



dt 



dtht-2n-2-\ 



.(n+D+i+^+D^+^y+^y 



* + ft- 



27rp 



n(,+ l) + ^-(2,+ l)f x +fe) 2 +(^)1 



t+n- 



pa-\-iz 
2irp 



(41) 



At this point it is convenient to introduce the complex quantity 

Si2=aS— x°i+i(r 2 — r 1 ) = — (pa+iz). (42) 

Then for brevity let 



2 12 « , •/ 

' = -2^ = 2^ + \ 



2?rp/ 



(43) 



Then l^oKl? since the absolute value, |zi 2 |, is the distance between 
the two poiuts, PJ and PI, in the reference sections, through which 
the helices hi and h 2 pass. The terms of the series in eq 41 may 
therefore be expanded before integrating. This gives 



dt 



+ L6 + 2^+(w + V2^/JJo WTo 

7i = l 

71 = 1 8=1 



71 = 1 
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+i;i»<»+'>2^-w]-^(^)} <«> 

*i — 1 O— 1 



n = l 



The first series in tliis equation is purely numerical and has the value 
.004, which is negligible, as are terms 10 times as large. 

The term containing the next series is found (by summing first 
with respect to n) to be 

-[-e+(#J+fe)>e»-«>- 

(71+ IV 

and write 

s=l s=l 

where 



and 
Hence 



= .154 . . . . , if 8=1. 

Hence introducing Biemann's zeta function 

CO 

r(«)»^for«>l, (45) 

»=1 

then the next to last series in eq 44 is 

8=2 

Similarly, to transform the last series in eq 44, we write 



8 = 1 
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where G n , i=B n , r and 

1 11 1 

Hence 

00 CD 

JV„ ,=2^^n=2f(2s-l), if «>1 

7» = 1 » = 1 

= .154, if 5=1. 
The last series in eq 44 is therefore 



Hence eq 44 becomes 



8=2 



il+i s - 1+ T ^1 27r ; 27r l_ 27r 2xp J 



+[-6-fi+(£j+(24) 2 ]- [1 ° g (1 -«- l0g *+** 



s=2 



where imaginary terms are to be discarded. 
This reduces to 

/ 1 +/ 2 =-l-g+ 1 r l (l- 1 J l )+1.85^„+.15^ 



s==2 *=2 

The definitions of I x and J 2 in eq 38 and 39 show that eq 37 is 

CO CO 

co a (X, a) 2 VVsos na l\^ T . n x , , /7 N1 cos na , r , r ,.,_. 
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Since we are considering the case where — 2^r<a<27r for which 
eq 46 applies, it is evident by reference to eq 28 and 28" that 



CO 

2 V Vos na _ 1 |a|/ jo: |\ 



(48) 



Also since X is small, as well as r 2 — r ly it is evident that k\ is small, 
so that by eq 27, neglecting infinitesimals 






cos na 



n- 



-f*i(*x) — 2 



a 
sln 2 



/ 2|sin§\ 



3^3 



. k\ II. «./ . ^^ 



2_1 

= 3 2 



sinf|(l-log|smf|)+[|-i|sinf|]log(g-|=) 



By eq 28' and 28" we may place 



so that 



. a 
sin ^ 



r\ 2%) 



-^.(W+^^^-f-^.-M)^-^^^] 



+ 



IHH0]*i&r « 9 > 



Using eq 46, 48, and 49, the eq 47 becomes 

«=2 1=2 

3a 2 ./ 2 V, z 
-4^- 2 ^/ l0g 2^ 

_H(' 1 _^r 1 _i og H(' 1 _y)4V^ ? 1i 1 4fe 



On replacing t by its value — ^~- and V r i r 2 by a, this gives 
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8 = 2 8 = 2 

The imaginary part is to be discarded. This equation is valid for 
|^|<C2*rp and for — 27r<«<27r. It shows that co a (A, a) contributes 
a correction term which is of the smallest order of magnitude retained. 
To express the first logarithmic term in eq 29 in terms of the variable 
2i2 defined in eq 42, we first let v=[r 2 — n|4-i(^2 — ^?)> so that, as in 
appendix 7 of paper 2 

_ Iog „_ .-?>_£_* i + 2ffi4(£)*, a M<2ip . 

s=l 

Since the real part of this is to be taken (this logarithm contains a 
real factor in eq 29), we may place 

log \v\= log \z 12 \ and @ {iv) 2s =^ (z l2 2 *,) 

since ^ (w) 2 '=^ [— (aS— a«) +i|r a — nl] 2 *, 

and this contains only even powers of — (a§— xj) and |r 2 — n). The 
eq 29 therefore becomes 



— TraJ-^ <l> 1 (k)+^-Jl—k 2 cos 2 1 

-(,FA C ° S 2; sm \*co8^-psin^log|_ jr j{» 

(51) 

where k t has been replaced by k and V r i r 2 by a in those terms which 
are of second order. In the term multiplied by 4:irN, we have replaced 
/ — u fi i ri— a+r 2 — a~] 
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6. SECOND-ORDER APPROXIMATION FOR m 

By use of eq 18, 21, 50, and 51 in the eq 15, we obtain the second- 
order approximation to m as the real part of eq 53 below. In this 
equation, we have placed 



\z l =x° 1 +i(r 1 —a 1 ) 
zi 2 =z 2 — zi, where \ 



(52) 



The complex quantities, Zi and z 2 , will be found to be the most 
appropriate ones with respect to integration over the sections of the 
wires. They are both referred to the centroid, Oi, of section Si as 
origin, but could equally well have been defined with reference to any 
other point, such as 2 . As shown by figure 1, the position of any 
point Pi(xi, 7*1) in section S t is represented by the vector or complex 
variable z X) and the position of any point P 2 (x 2 , ^2) in section S 2 by 
the vector or complex variable z 2 with reference to the same origin. 
Since P x and P 2 are in different wires, it is evident that \Z2-Zi\l2irp 
is less than 1. 

The equation for m takes the form 

m a (xl n; x° 2 , r 2 )=L s +2w x (l) 

where 

Ci=-\ "2.85+2^-2 log— 1 (53a) 

<7 2 =.85+27r<h-2 log ^ (53b) 

An 
<7 3 =.15 + 16^4+2 log— (53c) 

The constants A u A^ A iy and <fo are functions of k, which are defined 
in eq 19 and 8, respectively. 

-& -k-y l - ¥ cos2 1+ Kf) 2 ( cos i ) sin-1 (* cos l) 

<i-^M^ 4^^ \ (53d) 

where — 7r<a<7r. 
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The real part of the complex terms in eq 53 is to be taken, so that 
the principal branch of the multiple-valued function log (z 2 — z x ) may 
be understood, since this is indeterminate only to the extent of an 
additive constant which is a multiple of 2iri. With this understand- 
ing, this term and all the three other functions of (z 2 —Zi)/2tp are 
single-valued, analytic functions of either variable z x or z 2 in a region 
including the respective sections Si or S 2 . To define the other func- 
tions with general argument z (where z is later to be replaced by 
(z 2 —Zi)/2irp), we have the following series: 

m 

j[z) m log 2tt2- y~^ 2 *, when \z\<\, (53e) 

where log 2rz denotes the principal branch (which is real when z is 
positive). 

m 

fdz) - -2^^^ z 2 ', when |z|<l (53f) 

m 

fM * -^Y lTs-i ) g2 '"'' when l 2 l <L (53g) 



Obviously, 
so that 



«-2 

X(z)=2zf 2 (z), (53h) 

/ 1 (^)=2J^y 2 (20^ , =22/ 2 (^)--2J o y 2 (2 / )^ , . (53i) 



For integrating over sections not closely resembling circles, it will in 
general be more convenient to use these functions in the form of series. 
For circular sections, the use of the series for/ and / 2 may be avoided, 
as these are summable in terms of known functions as follows: 
Taking the logarithm of both sides of the equation 



sin tz 

7TZ 

gives 



HIO-I) 



n = l 

Since we are here concerned with the case |^|<1, we may place 



!og(i-|> 



._V!!!.JL, 

Si s n 2s ' 
so that 



, sin ttz ^ \\\\z 2 ' 1 _^ V^' V 1 

°^ irz ZjZj s n2s Zj s ' ^ 
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the inversion of the order of summation of the double series being 
valid because it is absolutely convergent. Since the zeta-function of 
Riemann is defined by 



this becomes 



^"S? ( when ) *>1, 



, sin tz VV(2s) 9 , , i i >- 

log -^- = -2j s Z > When M< L 



Comparing this with eq 53 e shows that if the logarithm denotes that 
branch which is real when its subject is positive, 

j(z)= log2sin7rz. (53 j) 

Similarly, we find that 

/ 2 (z)=2 7 z+log£j^, (53k) 

where Y=0.577216=Euler's constant. 

This may be found by taking the logarithm of the infinite product 
which represents the reciprocal of the gamma-function 



r(s+i) 

This gives 



--&** 



B =i\ nj 



When 



so that 



72 +lo g r (2+ i) = -2[log(i+|)-|} 

N<1 , _[ teg ( 1+ i)-|]_2(=^, 

*=*2 

yz+ log r(g+l) = V ( ~^' r(s) g ' 

9=2 

/ / 2s £i 2s— 1 

*=1 *«2 

1 , sin tz . 1 * , N 

=+| log r(i+^)+|iog r(i-«)+|/a(2)8 

since r(l + 2) r(l — z)=vzj&\tx tz. This proves eq 53k. 
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The principal term in eq 53 is the term Z„ representing the self 
inductance of a cylindrical current sheet, for as shown by eq 16, it 
contains the large factor N 2 or l 2 l{2irp) 2 . The correction terms of 
first order are two terms within the braces, the first being —2NJ on 
account of the large factor N. The second is the term containing 
A x with the large factor a/2irp. All other terms are of second order. 
The modulus, k y (complementary modulus k') is computed in terms of 
l=2wpN and D=2a, by k 2 =D 2 /(l 2 +D 2 ). The cylindrical radius a 
will be taken as the mean of the two mean cylindrical radii of the wires 
in the sense that a x and a 2 are the distances of the centroids of the 
reference sections from the axis. 

However, it is important to notice that eq 53 represents the mutual 
inductance, m a (zj r t ; x% fa), of any two helices, hi and h 2l with relative 
azimuth a, when the pair of coordinates (a§, r x ) and (x° 2 , r 2 ), which 
serve to specify hi and h 2 are such that \x 2 — x Q l \<C2irp i |aj<V a *id 
jr 2 — rx\ is of the order of magnitude of 2wp. More precisely \z 2 — Zi\<C2irp. 
It, of course, contains no reference to wire sections. Hence the 
constant a is, within certain small limits, still at our disposal. The 
constants ai and a 2 occur only in the first-order terms of eq 53, and 
there they appear in the sum, a 1 J ra 2 , which is 2a. 

The only thing that has been assumed concerning a in obtaining 
the approximations in eq 53 is that r x — a and r 2 — a are of the order of 
magnitude of the pitch. It is conceivable that for example the mean 
radius a x of the first wire section might be measurable with greater 
precision than a 2 . It might then suit our purpose to take a=a u which 
would then be used in determining the modulus k upon which all the 
constants i s , 2w z (l), C a , Ci, etc. depend. 

On placing a=0 in eq 53, it is found to be equivalent to the equa- 
tion which was used in paper 3 for the self inductance of a helical 
wire (neglecting D/4, which, as indicated above, is in general less 
than 1 part in a million) . 

IV. MUTUAL INDUCTANCE, M a , OF THE HELICAL 

WIRES 

In evaluating the integral of eq 4, which gives M a , there is an 
enormous saving of labor in adhering to the compex form until the 
integrations have been performed. Accordingly, we introduce the 
complex constant, z , which represents the vector from the centroid 
Oi of the reference section Sx to the centroid 2 of reference section S 2 . 
Its axial component is seen to be xl 2 of figure 1, which may now be 
designated by x . Its radial component is a 2 — ax. 
Hence 

Zo=Zo-\-i(a 2 —ax)=qe i0 , where — 7r<j8<7r, (54a) 

so that g=\z \ = -y/xl+(a 2 — a x ) 2 , and /3=tan -1 [(a 2 — aOMJ- The points 
of integration, P x and P 2 , are fixed respectively by the vectors z x and 
z 2 , both drawn from Ox as origin. The eq 53 then give 



x° 2 +i(r 2 -ax) = z +p 2 e%\- (54b) 
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Since each wire carries unit current, as implied by eq 3, it is evident 
that the constant terms in eq 53, namely, L s +2w x (l) + 2iraC ai when 
used in the integral 4 are merely multiplied by unity, whatever the 
shape of the sections Si and S 2 . 

All of the terms of eq 53, which are of second order, are then to be 
averaged first over one section, and this result then averaged over 
the other section. Since the points Oi and 2 are centroids of the 
reference sections, it is evident that the average value of x? over St 
is zero, and the average value of x 2 over S 2 is x . The average of 
(xi) 2 over Si may be denoted by k xx , and that of (x 2 — x ) 2 over S 2 by 
k' xz . Also the average of r t — k\ over S% is zero, as is the average of 
r 2 —a 2 over S 2 . We may denote by k vy the average of (n— a t ) 2 over Si, 
and by k' yy the average of (r 2 —a 2 ) 2 over S 2 . These averages multi- 
plied by the area of the section represent its moment of inertia. For 
the products of inertia, we use the similar notation, where k xv denotes 
the average over S x of £?(ri— a x ) and k f xv the average over S 2 of (x 2 — x Q ) 
(r 2 —a 2 ). 

Then k XXy k VV) k xy , etc. are all of the order of magnitude of (2wp) 2 or x%. 

The terms of first order in eq 53 must be treated more carefully. 
They are the two terms in 53 affected by the factors Aia/2wp and — 2N, 
respectively. For these we use the expressions 5a, and for the formu- 
lation in the case of general shapes of section we prefer to use rectangu- 
lar rather than polar coordinates in harmony with the definitions 
adopted for k xxi etc. 

The formula for the mutual inductance, M a , of two helical wires 
with general shape of reference section may then be written. 

M a =L s +2a> x (l) 

i i r> a I c x *ic xy ~r c Xt> Kx v ~\~c T Ac yy -\- C Tt K vv I 

+ 1<M 4. (5^F J 

where/, etc. denotes the double average taken in succession over the 
one section and then over the other, of the function of (z 2 —Zi)/2irp. 

The principal correction term comes from —2Nj, since this is the 
only term of first order in eq 55. 

In the case of circular reference sections k xy =0 and k xx —k vy =dl/16 
so this becomes 

+(^X*sf s I 86+2 **- sto '(i)] + 
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+i(?-0[ (1+2 ^ 
~- 2N f jit J !^K^0)[ cxiPi cos * i+ ( cri+ 2^) pi sin * 

+c X2 p 2 cos <j> 2 +( Crt+^j; ) p2 sin * 2 Jf ^ 

nd 2 
where dSi= p\dpid<j>i and Si=—r^' 

To perform the remaining integrations, consider the line average 
taken around any circular arc p x = const., of any monogenic function 
of 2i, multiplied by e in \. Considering n an integer, positive, negative, 
or zero, 



?\ 






yn Zi n ~ x F(zi)dz u since Zi = pie i<in 

This is a contour integral around a circle in the ^-plane, which 
encloses the point 2i=0. Hence by Cauchy's theorem 



*1* 



_ I e in *F(zi)d<l>i = 0, if n is a positive integer 

2 V-x 

=JF(0),if n=Q. 
and 

— I € ™*F{zdd4>i= — — - f ; if n is a positive integer, 

where 



Hence 



\ I F(z l )d<l> 1 =2F(0) (57a) 



Snow] 



I si 

V- 



cos n<j)iF(z l )d<t> i —- 
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Pl n F (n) (0) 



71/ 



sin 7Hl>iF(2i)d<t>i'> 



i Pl n F w (0) 



^ if 71 is a positive integer 
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(57b) 

(57c) 



Multiplying 57a by pidp u and integrating from pi=0 to pi^di/2, 
gives, on dividing by dJ/4, 



(5a) 



(58b) 

, if m is a positive integer. 

(58c) 



Similarly, from 57b and 57c, we find 

J J -^ „ cos n0^( 2l )^=-p3^ — ,- 

Taking F{z ] )=j[~ 2 ^—~ \ it is evident that the first average over the 
circle S% is 

J \ 27rp / J \2tvp) J \ 2irp / 

Letting z 2 denote p 2 e i(f>3 , the same theorem applies for averaging over 
the second circle, S 2 . 

Consequently, the four functions in eq 56, which are twice averaged, 
give 

Similarly, the formulas 58b and 58c suffice for evaluating the last 
integral in eq 56. 

The eq 56 then becomes 



229032—40 3 
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(59) 
where the real part is to be taken and 



z =Xo+i(a 2 — 0i) = 
Referring to eq 53, we find 



--qe* 



H^)= mlog ( 2 



-)=^log2sin[^±i|=^] 



2() 

2p 



fey- 



cot #°= 



I log 2 [cosh (2^)- cos |] (59a) 
.[sin|-isinh(^)] 



^/i( 



2 ? cosh/^^- 1 ) 



a-V-2f- ; (2s - 1 Y -i-'l : 







2x25/ J^ s \2irp/ 

f(2s-l)/ 



Xn 
COS — 



cos 2s/3 



(59b) 



(59c) 



<^--*g^(*F- <*•-»*■ 



2t^o 



+ lo£ 



■( 



r(i+ 



Zp \ 
2*ff/ 



r i- 



~2wpJ 



(59d) 



Using these real values in eq 59 gives the general formula for M a when 
the wires are circular. For practical purposes, this may be simplified, 
since any single-layer, bifilar winding which is to be used as an absolute 
standard of mutual inductance could be given the following char- 
acteristics. 

A. The two wires have practically the same diameter and the same 
treatment, so that the terms of eq 59 which contain the factor f 
will cancel each other (c xl =c X2 and c Tl =c T2 ). 

B. The mean radii a x and a 2 are approximately equal, so that we 
may place z =x in all except the principal correction term of first 
order, which is —2Nf (z /2irp). 

C. The helical windings are nearly equally spaced, so that z =x =7rp 
hVall terms except the principal correction term. 

In this case, by eq 59d, 
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<^MIF+^^ ; 



= 7— log 2 = — .116 



log^-==- log 2 = -.693, and from=53/ 

, (^> / <o=- 2 § f -^=-[± 8 i + ±C-^=- 1 )] 

The eq 59 then becomes 
M a =i,+2co I (Z) + 7rZ)((7 a +.12 + 1.31--.693^ 

+G+? )("*■- * | )- ivi °." <•** 2 f 2 '- cos l) 

+0-)('+^)(i)}< 6 °> 

It may be noticed that the term l-\-Dc r is +1 or — 1, according as the 
current distribution is uniform or "natural." The approximate 
value of the principal correction term within the brace is — N log 
4 = — 1.38629 N. 

When the initial sections of both wires are in the same axial plane, 
a=Qj the special case of eq 60 then is 

M.=L,+2»,ffl-,©{-.12+ilog(^)-^+g(l-^sfa-'f ; ) 

-0-O( i+ K>(^)' +wlo K cosh ^- cos f)l' <6Oa) 

When a= — ir, the special case of eq 60 is 
M„.-i,+2„,ffl-, D ji 1 o g (|f) + f + 2 Vs i lo i <i±£) 

-Kf- (1+1 ?4&) ,+ tf l °° i™^- ~fi ))• ( 6ow 

In this equation the lead wires may be taken in the same axial plane 
but on opposite sides of the cylindrical form, as shown in figure 2. 
As a guide in computing, it may be useful to remember that terms 
smaller than D/4 have been neglected in the derivation of formula 60. 
This corresponds to an uncertainty of about 0.1 in the value of all the 
terms within the brace which has the factor irD. 
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For the sake of comparison, the self inductance, L, of either of the 
two equal helical wires (paper 3) may be put in the form 

Z=Z,+2o;,(Z) 

+ 7rZ>{2iv[log ^-.89473]+! log 

wfr k'/, k' 



D 

2p 



3 



W-T** 1 *)}* < 61 > 



where o) x (l) in this equation and in eq 55, 56, 60a and 60b is computed 
by eq 22. 

V. EFFECT OF LEAD WIRES 

When the lead wires are each parallel and close to a generator of 
the cylindrical form, it is simpler to use, in eq 60, instead of the term 
2(ti x (l), the term 2o) x (l, a), which takes account of all axial currents in 
lead wires and helices together. In paper 1 eq 42 to 46 are given the 
relations which in the case r 1 =r 2 =D/2 are equivalent to 

2«,a, «) =2 { v(?-0 +B(a) - -JP+B& 

-I 2 , „,„ DU+JP+B&] j . . 

7?w log 2+Uog R(«)uU*+mv m 

where 

R 2 (a)=r l 2 1 -2r h r l2 cos a+rf 2 (63) 

the distance of the lead wires from the axis of the cylindrical form 
being r h and r* 2 , where r h — a and r H — a are small, so that when 
o:=±7r, R(±tt) is practically the same as the mean diameter, D, of 
the coils. 

The term 2zs x (l, a) is to be substituted for (not added to) the term 
2u x (l) in eq 60 and its special cases. 

VI. MODIFICATION REQUIRED WHEN THE REFERENCE 
SECTIONS ARE NOT QUITE CIRCULAR 

In this case, we first find the position of the centroids of the sections, 
thus determining a u a 2y and x 0) so that the points 0i and 2 are deter- 
mined. With these points as centers, we describe two circles whose 
diameters, d x and d 2 are determined by 

*(f V=Si and *(^) 2 ==S 2 . (64a) 

The equations of the boundary curves of Si and S 2 may then be put 
in the form 

Pi=~[l+ei(d>i)] and p 2 =^ 2 [l + e 2 (<fe)], (64b) 

where e x and e 2 are functions of #1 and fa, which must take on positive 
and negative values. We assume that the order of magnitude of 
|ei| and |€ 2 | is not greater than^that of p/a. 
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We may at first consider S 2 to be exactly circular. The section Si 
differs from the area of a circle of diameter d\ by the area of a narrow 
strip whose radial breadth is 

Since the centroid of circle is at U and also the centroid of 8% is at 0\ } 
it follows that the centroid of this narrow strip is at 0\. Hence we 
have imposed four conditions upon the small function €i(0i). 

(d Y 
i? ) ei(<l>i)d<l>i, these four 

conditions are 

\*i(4>i)\^C-> (where G is of the order of magnitude of unity) (65a) 

f*«i(*iy*i = (65b) 

€i(&)cos0id0i=O (65c) 

J €i (&) sin fafox = (65d) 

The eq 65b states that d\ has been so chosen that the area of the 
noncircular section, Si, is vl %r J 

Equations 65c and 65d state that O x is the centroid of the section S% 
and also of the circle. 

The function e x must be a periodic function with period 2ir. Hence 
it is developable in a Fourier's series in the form 

CO 

€i(0i)= V(<7* cos nfa+Dn sin ufa), (66) 

11=2 

where the term n=l must be absent because of eq 65c and 65d. The 
constant term is absent because of condition 65b. 

As explained in the derivation of eq 59, the only effect of this 
deviation of Si from exact circularity will consist in a correction (of 
second order) that must be made to the principal correction term 
of eq 59, which is the term —2NJ(z /2tp) within the brace. 

Instead of the termj(z /2Trp) in eq 59, we must place 

r rdSt r rds 2 jzo-zA c rds 1 J*>-zi\ 
J J TS J J T*f\-2^rr) J srA"2iF/ 

where the integral is now taken over the section Si, whose boundary 
curve is given by eq 64b. When this same integral was taken over 
the circular section, its value was found to be j(z /2irp). Hence the 
correction term to be added to eq 59 is equal to 1/Si times the integral 
oiJ(z Q — Zi)/2irp) taken over the area of the narrow strip. This gives, 
since €i is small, the correction 
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where the integral is taken around the circle px—di/2. When A x is 
found, the term —2Nf(z /2irp) in eq 59 must be replaced by 

— 2iV Ai+yfo - ^") * This will make a change in the formula 59 

by an amount which cannot be larger than a second-order term. 

If, in eq 67, we insert the series 66 and make use of eq 57b and 57c, 
we find 

^yu^m^y^y (68a) 

w=2 

Similarly, 

n=2 

when the boundary curve of S 2 is of the type of eq 64b and 

00 

€2(<fe)= 2 (<% cos n<t> 2 +D' n sin nfa). 

w=2 

Therefore, the correction is made to eq 59, 60, 60a, and 60b by adding 
within the brace which is multiplied by 7rZ>, the term — 2iV r (A 1 +A 2 ). 
Perhaps the case of most practical interest is the elliptic section 

!+^f^=l, (69) 

or in polar coordinates 

A 1 B l 



Pi= 



If (Ai—B 1 )/A l is small, this becomes, neglecting (A—Bx) 2 /A 2 , 

Pl = V3A [i+^gf 1 cos 2*i] • 

The diameter, d h of a circle whose area is equal to the area tAiBi of 
the ellipse is given by di/ 2 =^AiB , hence the boundary curve is 

Pi=f [l+^f^ cos 2*] 1 (70a) 

so that, on reference to eq 64b, we find 

^ifoiH^p^ 1 cos 2*. (70b) 

The only Fourier's coefficient in eq 66 which is not zero is C% % which is 
equal to (Ai—B^/di. Hence, by eq 68a, 
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From eq 53 j, we find ... . f*\&)=*- 



2 ^0 

sin 2 — 



2p 



Taking the real part 

,(a 2 —a l \ x 



*■&>-- C o Sh (^) + c <|)' 



so that 



a _ _ 1/4.-1^ Y C0Sh V"^rh C0S 2i > 

2^ 2p 

The correction, A 2 , by change of subscripts, takes account of the fact 
that the reference section of the second wire is an ellipse, as in eq 69, 
with semiaxes A 2 and B 2 . In these equations A\ may be greater or 
less than B\. Hence with both wires of elliptic sections the correction 
to be added to eq 59 within the brace is the real quantity 

-2N(A 1 +A 2 ) 

7 a<i — Qj\ Xq 

2p 2p 

Since, by hypothesis, the order of magnitude of (A 1 —B 1 )/d u etc., is 
not greater than that of p/a, it is evident that the correction term 71 
is not greater than second order. If (A—iB^/di is of the order of 
magnitude of p 2 /a 2 , the correction for ellipticity is entirely negligible 
in the formulas 59 to 60b. 

In the eq 60 and its special cases 60a and 60b we have placed 
Xo=irp and a 2 — a 1 =0 in the second-order terms of eq 59. Hence for the 
eq 60, 60a, and 60b, the correction term 71 reduces to 

-*»<^-mi^xi)wxi)*j »« 

We have idealized the actual helical wire, w% } as one which is cut 
in a circular section by a plane through the axis of the cylindrical 
form. This is the circle shown in the reference plane, or xy plane, 
of figure 1. However, a more natural view of the wire is one which 
is cut in a circular section by a plane normal to its central filament. 
The reference section of such a wire would be slightly elliptical. If 
its semiaxes are A\ and B u it will be shown that the ratio (Ai-—B 1 )/di 
is equal to p 2 /2a 2 , so that by eq 71 the correction is of third order and 
is therefore not to be made to the formula here given, since we have 
neglected an unknown number of terms of third order. 
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To examine this second view of the helical wire, consider the three 
eq 1, which may be written (dropping subscripts) 

X=x°+p I I — 2irp<Cx°<i2Trp} 

y=r cos 0>, where j— co<^0< oo I (72) 

z—r sin J [ 0<r<oo J 

Since the assignment of the values of x°, r, and $ uniquely determines 
the position of a point in space, they constitute a system of space 
coordinates which are not orthogonal. With this specification of their 
ranges, they are uniquely determined when x, y, and z are given. 
Any equation of the form £°=0(r, 0) represents a surface. In partic- 
ular, the equation independent of 

z°=<Kr) (73) 

represents the surface of a helicoidal wire. This surface cuts the 
cylinder r= constant in a helix whose equations would be given by 
eq 72 on replacing x° by 0(f), and holding r constant, while is varied. 

The curve in which the surface of eq 73 is cut by a plane, 0=c= con- 
stant, is congruous with the curve in which the surface is cut by any 
other plane, 0=Ci= constant. The two curves may be made to coin- 
cide by giving one a translation in the x direction, followed by a 
rotation about the x axis. 

Taking as a representative cutting plane, the plane 0=0 gives a 
curve whose rectangular equations are by eq 72 

z=4>(y) (74a) 

2=0. (74b) 

Hence if we have the equation of the surface of any helical wire in 
the form of eq 73, we may find the curve bounding its reference section 
of figure 1 by replacing x° by x and r by y in eq 73. 

Next, suppose the equations of any closed space curve are given 
in rectangular coordinates in the form 

f l (x f y ) z)=0 (75a) 

f 2 (x,y,z)=0. (75b) 

Let it be required to find the equation of the surface of a helical wire 
which passes through this assigned curve. 

Writing the equations of the given curve in helicoidal coordinates 
£°, r, and 6, they are, by eq 72, 

r l (x°+p0, r cos 0, r sin 0) = (76a) 

f2(3?+p0, r cos 0, r sin 0) = 0. (76b) 

Eliminating between these two equations, we obtain a relation 
between x° and r, which, if solved for x°, takes the form of eq 73. It 
is therefore the required equation of the surface of the helicoidal wire 
but expressed in helicoidal coordinates. Its reference section is 
found by placing x°=x, r=y, as in eq 74. 

Such elimination of is seldom possible, so that the reference 
section may be considered as a curve defined in parametric form 
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by placing x°=x and r=y in the two eq 76a and 76b. The variable 
line parameter is 0. 

Consider next another rectangular system of coordinates obtained 
by rotating the x, y, z system of figure 1 about the y axis, through a 
positive angle, y. Any point, P, in space has the rectangular coor- 
dinates x, y, z, in the old system and x' , y', z' , in the new, where 

x'=x cos y—z sin 7 

y'=y 

z'=x sin 7+2 cos 7. (77) 

Suppose that we desire the equation of the surface of a helical wire 
whose section by the plane z' = Q is a circle of diameter d with center 
at x' = 0, y'=a, z' = Q. The equations of this circle (as a space 
curve) are 

x*+<3j'-a)*=*j (78a) 

s'=0 (78b) 

or by eq 77 

d 2 
(x cos 7—2 sin y) 2J r{y—a) 2 =-r (79a) 

x sin 7+2 cos 7 = 0. (79b) 

The two latter correspond to eq 75a and 75b. Inhelicoidal coordinates, 

these become 

d 2 
[(x°+p$)oos 7— r sin sin y] 2 +(r cos d—a) 2 —j=0 (80a) 

(x°+pB) sin y+r sin cos 7 = 0, (80b) 

which are special cases of the general eq 76a and 76b. The relation 
between a; and r resulting from the elimination of between 80a and 
80b would be the equation of the surface of a helicoidal wire which 
is cut by the plane z f = (i. e. z— — x tau 7) in the circular section de- 
fined by 78a and 78b. Instead of attempting to eliminate 0, it is 
better to retain it as a parameter. The equations of the boundary 
curve of the reference section of this wire is found, as explained above, 
by replacing x° by x and r by y. They are 

d 2 
[(x+pd) cos y—y sin sin y] 2 -\- (y cos d—a) 2 —-j = (81a) 

(x+p6) sin y+y sin cos 7=0 (81b) 

(together with z=0). 

The second view of the wire corresponds to 

a T) 

cos 7= . and sin 7== / > , J (82) 

iar+p 2 -ylar+V 
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which, as shown by 72, makes the plane z'=Q normal to the central 
filament of the wire at z=0, y=a, 2=0. Using eq 82, eq 81 may be 
written 

. /asin0\ , N 

x = -pe-y^—— J (83a) 

{a 2 sin 2 6+p 2 )y 2 -2p 2 a cos e.y+p 2 (a 2 -~)=0, (83b) 

which are the parametric equations of the boundary curve of the 
reference section of the wire. The last was obtained by eliminating 
x between 81a and 81b. On solving this equation for y, we obtain 
the following parametric equations for the boundary of the reference 
section of this wire: 

»~!"-< £ f L ') <84«> 



y= 



/ a sin fl V 
\ P / 



(84b) 



where p=d/2. 
From this we find 



dy__ ay sin d 
dx~~ p(a cos d— y) 



The section is a closed curve symmetrical with respect to the line 
x=0. 

The real values of y correspond to values of the parameter 6 in the 
range — o <0<0o, 



where 



o =sin"" 1 



PP 



a a 2 +p 2 —p 2 



If we now introduce the facts that p/a and p/a(=d/2a) are small, 
it is found that 



*-£['+*-($)] 



Hence in eqs 83a and 83b we may expand sin 6 and cos 6 in powers of d 
This gives, neglecting terms smaller than p 2 ja 2 , 



?(■+&» 






On eliminating 6 between these equations, we find for the boundary 
of the reference section, the equation 
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where 
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x 2 (y-a) 2 _ 






A 2 ' B 2 l$ 




A= 


A-B f 

d 2a 2 ' 


d 

= 2' 



so that 



as stated above. 

Since terms smaller than second order have been neglected, the formula 
here obtained for the mutual inductance of the two helical wires has 
the same validity for either of the two modes of defining the helical 
wire. This applies to formulas for mutual inductance and selfinduct- 
ance obtained in the earlier papers quoted above. The conclusion 
is important because the absolute determinations of the ohm and of 
the ampere which have been made at this Bureau during the last 
decade are based upon the use of these formulas, in which the terms of 
second order have not been negligible. However, if the wire is found 
to have an elliptic section whose elhpticity is a small quantity of the 
first order, then the correction 71 must be applied. 

VII. SUMMARY 

The formula for the mutual inductance of the two helical wires of 
wire diameter d, which are wound bifilarly in a single layer on a 
circular cylindrical form is given by eq 60. In this formula the 
azimuthal angle, a between the initial ends of wires w x and w 2 lies 
between — t and ir, inclusive. When a=— w, the central filament of 
wire ^begins at 2 in figure 2; when a=0, it begins at 2 ; and when 
a= + ir, at 0' 2 '; the central filament of wire w x beginning at 0i in all 
cases. The total number of turns in each helix is N. The cylindrical 
radii, a x and a 2 , of the central filaments are considered so nearly equal 
that a 2 —a 1 /27rp is small where 2irp is the pitch of the windings, the 
axial length, I, being 2-n-pN. The principal term of eq 60 is i„ which 
represents the self inductance of a circular cylindrical current sheet 
of length I and radius a=(a 1 +a 2 )/2=Dl2. This sheet has angular 
component of current only. 

The term 2ca x (l) in eq 60 represents the effect of axial components 
of current in the helices. It is a relatively small term of second order 
which may be computed by eq 22, but which should be replaced by 
27o x (l) given by eq 62 when the lead wires are close to a generator of 
the cylinder, as this takes into account the effect of axial currents in 
lead wires and helices together. 

The principal correction term in eq 60 contains the large factor, N. 
In this term, x is the (positive) axial component of the line from 0i to 
Q2 in figure 1. In all other correction terms, x has been placed equal 
to wp and a 2 —a x equal to zero, but allowance is made in this term 
for the possibility that measurements will show slight deviations 
from these values. The value x =wp corresponds to the fact that 
the windings of one wire are equally spaced between those of the other. 

In section VI, a modification of the formula is given to apply when 
the reference sections of the wires are only approximately circular. 
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VIII. LIST OF SYMBOLS 

0i= mean radius of primary wire, W\. 
a 2 =mean radius of secondary wire, w 2 . 
d t = diameter of wire W\. 
d 2 = diameter of wire w 2 . 
27rp=the pitch of each winding. 
N= number of turns in each winding. 
l=2TrpN=the axial length of each winding. 

x =7% of figure 1 = axial distance from center of wire Wi to center of 
wire w 2l both centers lying^ in a plane through the x axis. 
When windings of one are midway between those of the other, 
Xo=irp. It is assumed that x is near irp. q and defined 
by 2 ==z +i(a 2 — ai) = £e lV? , where — 7r/2</3<7r/2. 
a=the difference in azimuth of the initial points of wires w x and 
w 2 , both wires having the same end planes when a=—ir. It 
is assumed that |a|<y» # 
c x , Cr , Cx t c T = partial derivatives of the logarithm of the magnitude 
of the current density, for points on the respective central filaments 
(eq 5b). For uniform current distribution they are zero; for the 
natural distribution 
c x =c x =0 and c Ti = — l/a h c r —\\a 2 . 
D=2a=a 1 +a 2 . 

D 2 2 

Jc 2 = p>jj2 an d V =1— k 2 . 

K and E are the complete elliptic integrals with modulus k. 

i,=the self inductance of a current sheet of length I and diameter 
D (eq 16). 

M=the mutual inductance of the helical wires, w x and w 2 . 

00 

f(2s-l)=]>\^i(= 45). 

71 = 1 

w<l>^l[K-E]-kK(=S). 

C a defined in eq 53 d. 
Conditions A, B, and C preceding eq 60 give the assumptions by 
which eq 60 becomes a special case of eq 59. 

The term 2w x (Z) in eq 59, 60, 60a, and 60b, may be computed by 
eq 22 and 23. It represents the effect of the axial components of 
current in the two helical wires, w% and w 2 . The term 2co r (7, a) may 
be computed by eq 62 and 63. When the lead wires are close to the 
cylindrical form, as shown in figure 2, the term 2c3 x (Z, a) should 
replace 2u x (l) in eq 59 to 60 to take account of all axial components 
of current in helices and lead wires together. 

The "helical wire" is assumed to be built up of helical filaments 
(all with the same pitch) in such a manner that any plane through the 
x axis cuts the wire in circular sections. 
i 

Washington, December 19, 1939. 



